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Abstract 

The single-particle spectrum obtained from the relativistic mean field (RMF) theory is used to 
extract the shell correction energy with the Strutinsky method. Considering the delicate balance 
between the plateau condition in the Strutinsky smoothing procedure and the convergence for the 
total binding energy, the proper space sizes used in solving the RMF equations are investigated in 
detail by taking 208 Pb as an example. With the proper space sizes, almost the same shell correction 
energies are obtained by solving the RMF equations either on basis space or in coordinate space. 
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The liquid drop model of the nucleus was historically the first model to be proposed as 
an explanation of the properties of the nucleus. Then the occurrence of the magic numbers 
has been one of the strongest evidences for the nuclear shells. The shell effect presents 
a fluctuation in the binding energy, the so-called shell correction energy, which can be 
supplemented to the liquid drop model to improve the descriptions of the nuclear masses. 
One of the most practical and effective methods for extracting the shell correction energy 
was proposed by Strutinsky l|, |2| . The Strutinsky method has been widely applied to the 
calculations of the masses, shapes and fission barriers of nuclei {3, 4]. 

Two key points in the Strutinsky smoothing procedure are the so-called plateau condition 
as well as the reliable single-particle spectrum. 

On one hand, the plateau condition means that the shell correction energy should be 
insensitive to the smoothing range 7 and the order of generalized Laguerre polynomial M , 
since these two variables are only parameters in the smoothing procedure which have no 
physical meanings . The plateau condition has been investigated in the past several 
decades. The Strutinsky method was first combined with the Nilsson model. In this case 

n 

a good plateau could appear due to the infinity of the potential pj. However, no plateau 
appears in the realistic finite potential because of the discontinuity of the level density at 
the upper edge of the potential A possible solution for this problem is to complete the 
spectrum outside the bound states region, using the resonances in the continuum region 6| 
or using the positive eigenvalues obtained by diagonalizing the Hamiltonian in a harmonic- 
oscillator basis {4], Q]. The first method is cumbersome, especially for deformed nuclei, 
although the plateau could appear if resonances up to 60 MeV are included {§]. The second 
one is more practical, however, the shell correction energy and even the appearance of 
the plateau condition depend on the number of harmonic oscillator shells included in the 
basisP, ~\. This requires a delicate adjustment to the shell number N . In the non-relativistic 
framework, the plateau condition for 208 Pb has been checked with the phenomenological 
folded Yukawa potential^, 4], where the shell correction energy was found to be stable in 
the range of 8 < Nq < 13, thus Nq = 12 was recommended and a good plateau was found. 

On the other hand, the self-consistent relativistic mean field (RMF) theory 
received wide attention due to its successful description of lots of nuclear phenomena[9| 



which has 



111 ], could provide the single-particle spectrum microscopically. With this single-particle 



spectrum thus obtained, the shell correction energies in superheavy nucleus have been inves- 
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tigated and it is found that the shell correction energy at the saddle point is too important 



to be neglected 



121 ]. which will influence the synthesis of the superheavy nuclei 
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The RMF equations can be solved by diagonalizing the Hamiltonian in a harmonic- 



oscillator basis space 



151 ] or by solving the coupled channel differential equations within a 



finite box [0, R] in coordinate space la], meanwhile the single-particle spectrum outside the 
bound states region can also be obtained. In this way, the level density in the continuum 
depends on the space size, i.e. the shell number for fermions Nf (bosons Nj,) or the box 
size R. Therefore, the space size used in solving the RMF equations is one of the crucial 
quantities for the Strutinsky smoothing procedure. Meanwhile, the choice of the space size 
should guarantee the convergence for the total binding energy. 

In this Letter, based on the single-particle spectrum obtained from the RMF theory, 
the shell correction energy will be extracted by the Strutinsky method. The efforts will be 
focused on determining the proper space size used in solving the RMF equations via the 
delicate balance between the plateau condition and the convergence for the total binding 
energy. 

The basic ansatz of the RMF theory is a Lagrangian density where nucleons are described 
as Dirac spinors that interact via the exchange of several mesons (the scalar cr, the vector 
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u, and isovector vector p) and the photon [8j. 

The classical variation principle gives equations of motion for the nucleon, mesons and 
the photon. As in many applications, we study the ground-state properties of nuclei with 
time reversal symmetry, thus the nucleon spinors are the eigenvectors of the stationary Dirac 
equation 

[ a . p + V (r) + /3{M + S{r))\ ^ = e^, (1) 
and equations of motion for mesons and photon are 

(-A + ml) a(r) = -gaps(r) - g 2 v 2 - gs^ 3 , 
(-A + ml)u\r) = g w p v {r), 
(-A + m 2 p )p°(r) = g p p 3 (r), 
—AA°(r) = e Pp (r), 

where u° and A are timelike components of the vector u and the photon fields and p° the 
third component of the timelike component of the isovector vector p meson. Eq. ([I]) and 
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Eq. (|2j) are coupled to each other by the vector and scalar potentials 

S(r) = g a a, 

and various densities 



i=i 

A 

i=l 

Pair) = ^2i/j}T 3 i/ji, 

i=l 

A 1 _ 



(3) 



(4) 



P. 

i=l 

The above RMF equations can be solved in basis space or coordinate space as men- 
tioned above. The single-particle spectrum thus obtained will be used to calculate the shell 
correction energies with the Strutinsky method. 

The Strutinsky method is based on the assumption that the realistic shell correction 
energy can be extracted uniquely from a mean field model. The shell correction energy, 
representing the fluctuating part of the binding energy, is defined as the difference between 



the total single-particle energy and its smooth part 



5] 



N(Z) ~ x 

E shcll = E-E=J2d- 2 / C9(c)de, (5) 
i=i 

where N (Z) is the neutron (proton) number, is the single-particle energy. The smoothed 

f X 

Fermi level A is determined by the conservation of the particle number, N(Z) = 2 / g(e)de. 

»/ — oo 

The smoothed level density g(e) takes the form, 

1 _ 2 

with the smoothing range 7. The folding function is usually taken as f(x) = —^e x P(x), 
where P(x) is an generalized Laguerre polynomial L^ 2 (a; 2 ). In practice, a cutoff is needed 



in the summation of the single-particle energy, and the shell correction energy is stable when 
the cutoff energy is large enough. 

The Strutinsky method, as shown above, has two additional variables, the smoothing 
range 7 and the order of generalized Laguerre polynomial M. Since neither 7 nor M has 
physical meaning, the value of the shell correction should be insensitive to these quantities 
within a certain range of values, i.e. the so-called plateau condition, 

9E shc \\ _ dE shcn _ 

This forms an important criterion for the feasibility of the Strutinsky method. 

Taking the nucleus 208 Pb as an example, the neutron single-particle spectrum is calculated 
by the RMF theory with the parameter set PKljl^]. In the Strutinsky smoothing procedure, 
the unit of the smoothing range 7 is adopted as htv = 41A" 1 / 3 (1 + MeV, and the 

energy cutoff is 30 MeV. 

First of all, solving the RMF equations in basis space, it is found that the shell correc- 
tion energies are almost independent on the shell number for bosons Nb, which could be 
understood as follows. The shell correction energies are determined by the single-particle 
level density. Since the single-particle energies are the eigenvalues of the Hamiltonian matrix 
for fermions (nucleons), changing the shell number for fermions Nf will change the single- 
particle level density in the positive energy region, whereas changing Nb doesn't influence 
it. Thus, the shell correction energies are not sensitive to iV&. In the following investigation, 
N b = 20 is fixed as a proper number [3] , and the efforts will be focused on the influence of 
Nf on the shell correction energies. 

In the upper panel of Fig. [T], the neutron shell correction energies as a function of the 
smoothing range 7 for the nucleus 208 Pb are shown, where the RMF equations are solved 
in basis space with different shell numbers Nf = 10 ~ 20 for fermions, and with fixed shell 
number Nf, = 20 for bosons. A well-pronounced plateau is seen in the case of Nf = 12. It 
is also found that the plateau vanishes when Nf > 12. The reason is similar as the case 
not including the positive energy part [6(. When Nf increases, the single-particle levels 
in the positive energy region become much denser than those in the bound region. The 
discontinuity in the single-particle level density makes the plateau vanish. 

Another important criterion for choosing the shell numbers is that Nf should be large 
enough to achieve the convergence for the total binding energy. According to this point, 
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Nf > 12 was proposed [181 ] . Specifically, in the present calculation, the accuracy for the 
total binding energy is 0.04% when Nf = 12 is chosen. 

Taking the balance between the above criteria, the size of basis space used in solving 
the RMF equations should be taken as Nf = 12 and N^ = 20, which is in accordance with 
the space size for folded- Yukawa potential 4j. The neutron shell correction energy reads 
-11.92 ~ -11.74 MeV within the range 1.2 < 7 < 1.8. 

In the lower panel of Fig. [TJ the neutron shell correction energies as a function of the 
smoothing range 7 for the nucleus 208 Pb are shown, where the RMF equations are solved 
in coordinate space with different box sizes .Rbox = 11 ~ 19 fm. The plateau condition is 
well fulfilled in the cases of i?box < 15 fm. Especially, the shell correction energies are stable 
between -Rbox = 11 fm and -Rbox = 13 fm. It is also found that the plateau vanishes gradually 
as -Rbox > 15 fm. The reason is the same as that in the basis space. 

On the other hand, -Rbox should be large enough to achieve the convergence for the total 
binding energy. It is found that the accuracy for the total binding energy is 0.01% when 
-Rbox = 13 fm is chosen. 

Therefore, based on the above two criteria, the box size -Rbox = 13 fm is recommended for 
solving the RMF equations in coordinate space. The neutron shell correction energy reads 
-12.00 ~ -11.80 MeV within the range 1.2 < 7 < 1.8. 

In order to verify the choice of the above space sizes, the dependence of the plateau 
condition on the order of generalized Laguerre polynomial M will be investigated. 

In Fig. [21 the neutron shell correction energies as a function of the smoothing range 7 
for the nucleus 208 Pb are shown. The RMF equations are solved in basis space with the 
shell numbers Nf = 12, iVj, = 20 (the upper panel), as well as in coordinate space with the 
box size .Rbox = 13 fm (the lower panel). The four different curves correspond to the order 
M = 2,3,4,5 of generalized Laguerre polynomial, respectively. It is found that there is no 
clear plateau in the case of M = 2, whereas well-pronounced plateaus appear in the cases 
of M = 3,4,5, and they are almost identical. Thus, the shell correction energies calculated 
with the above selected space sizes are independent on the order M of generalized Laguerre 
polynomial within a reasonable range. Furthermore, M = 3 is an optimal value for the 
present calculations, since the optimal value of M is reached as soon as the shell correction 
energy remains constant when M is increased {5]. 

In summary, considering the delicate balance between the plateau condition in the Struti- 
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nsky smoothing procedure and the convergence for the total binding energy, the proper space 
sizes used in solving the RMF equations are investigated in detail by taking 208 Pb as an ex- 
ample. Solving the RMF equations in basis space, the shell number for fermions Nf = 12 
is recommended, while it is found that the shell correction energies are almost independent 
on the shell number for bosons N^. Meanwhile, -Rbox = 13 fm is recommended when the 
RMF equations are solved in coordinate space. The above choices of the space sizes lead to 
almost identical neutron shell correction energies i£ S heU ~ —11.90 MeV in both cases, and 
are verified by the independence of the shell correction energy on the order M of generalized 
Laguerre polynomial around M — 3. 
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FIG. 1: Neutron shell correction energies as a function of the smoothing range 7 for 208 Pb 
calculated by RMF theory with PK1 parameter. The RMF equations are solved in basis space 
with the shell numbers Nf = 10, 12, 14, 16, 18, 20 for fermions and the fixed shell number Nf, = 20 for 
bosons (the upper panel), as well as in coordinate space with the box sizes Rhox = H, 13, 15, 17, 19 
fm (the lower panel), respectively. 
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FIG. 2: Neutron shell correction energies as a function of the smoothing range 7 for 208 Pb calculated 
by RMF theory with PK1 parameter. The RMF equations are solved in basis space with the shell 
numbers Nf = 12, Nf, = 20 (the upper panel), as well as in coordinate space with the box size 
-Rbox = 13 fm (the lower panel). The four different curves correspond to the order M = 2, 3, 4, 5 of 
generalized Laguerre polynomial, respectively. 
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